We study the expansion of the ionization and dissoci- 
ation fronts (DFs) in a radially stratified molecular cloud, 
whose density distribution is represented as n(r) oc r~ w . 
We focus on cases with w < 1.5, when the ionization front 
is "trapped" in the cloud and expands with the preceding- 
shock front. The simultaneous evolution of the outer pho- 
todissociation region (PDR) is examined in detail. First, 
we analytically probe the time evolution of the column 
densities of the shell and envelope outside the HII re- 
gion, which are key physical quantities for the shielding of 
dissociating photons. Next, we perform numerical calcu- 
lations, and study how the thermal/chemical structure of 
the outer PDR changes with different density gradients. 

■ We apply our numerical model to the Galactic HII re- 
gion, Sharpless 219 (Sh219). The time evolution of the 

■ column densities of the shell and outer envelope depends 
\ on w, and qualitatively changes across w = 1. In the cloud 

■ with w < 1, the shell column density increases as the HII 
, region expands. The DFs are finally trapped in the shell, 

and the molecular gas gradually accumulates in the shell. 
| The molecular shell and envelope surround the H II re- 

■ gion. With w > 1, on the other hand, the shell column 
| density initially increases, but finally decreases. The col- 

■ umn density of the outer envelope also quickly decreases 
| as the HII region swells up. It becomes easier and easier 

for the dissociating photons to penetrate the shell and en- 
velope. The PDR broadly extends around the trapped H II 
region. A model with w = 1.5 successfully explains the ob- 
\Q | servational properties of Sh219. Our model suggests that 

■ a density-bounded PDR surrounds the photon-bounded 
^ I HII region in Sh219. 

Oh; 

' | Key words. Stars: early-type - Stars: formation - ISM: 
^ ■ HII regions - ISM: kinematics and dynamics - ISM: 
' molecules 
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1. Introduction 

H II regions are one of the basic elements of the interstellar 
medium and have been studied by many authors. Roughly 
speaking, the theoretical modeling has been advanced 
from two complementary standpoints. Some studies have 
modeled the detailed ther mal and chemical st ructure of 
the static HII region (e.g.. iFerland et aDll998h and pho- 
todissociation region fPDR. e.g.. iTielens fc Hollenbachl 
1198.4 lAbel et aljbTiol . We can compare these models 
with the observed line-strengths and their ratios to probe 
the physical state of the target region. Other studies have 
focused on the d ynamical tim e evolution of the expanding 
HII region (e.g.. lYorkelll98qh The density and tempera- 
ture actually evolve as the HII region expands, and the 
gas flows across the ionization front (IF). A shock front 
(SF) sometimes precedes the IF owing to high pressure 
in the HII region. The dynami cal evolution is diffe rent in 
different density distributions ijFranco et al.lll990l here- 
after FTB90). The ionized gas sometimes blows out from 
the m olecular cloud as a "champagne flow" l|Tenorio-Taglel 
Il979h . 

Confronting some recent observations, however, we re- 
quire up-to-date modeling which includes both the fluid 
dynamics and c onsistent therm al/chemical structure (e.g., 
IFerland! 120031 iHennevI l2006|) . For example, the "col- 
lect and collapse" by the expanding HII region is one 
of th e dynamical triggering processes of star formation 
(e.g.. lElmegreen fc Ladall977tlElmegreerJl998il . Recently, 
high-resolution observations have revealed the detailed 
structure of the shell around some pc-scale H II regions 
foeharveng et aD l2003al iDeharveng. Zavagno k, CaplarJ 
120051 IZavagno et al.ll2006|) . The dense shell-like struc- 
tures are traced by the molecular and dust emission just 
around the ionized gas. Young clusters including OB stars 
are often embedded in fragments of the shell, which is 
along the "collect and collapse" scenario. In our previ- 
ous papers l|Hosokawa & Inutsukal l2005. 2006al hereafter, 
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Papers I and II), we have focused on the fact that a cold 
molecular rather than warm neutral shell surrounds the 
H II region. We have analyzed the dynamical expansion 
of the HII region, PDR, and the swept-up shell, solv- 
ing the \JY{hv > 13.6 eV)- and FVY{hv < 13.6 eV)- 
radiative transfer, thermal and chemical processes in a 
time-dependent hydrodynamics code. Our numerical cal- 
culations have shown excellent agreement with some ob- 
servational properties of a "collect and collapse" candi- 
date. The PDR is quickly trapped in the shell, and a cold 
and dense molecular shell forms around the H II region in 
a homogeneous ambient medium. 

In this paper, we examine how such basic time evolu- 
tion changes in an inhomogeneous ambient medium. Since 
massive stars form in a dense region in the molecular 
cloud, such a situation will be more realistic than expan- 
sion in the homogeneous medium. Actually, observational 
features of some H II regions are very different from those 
of the "collect and collapse" candidates. For example, the 
Galactic HII region, Sharpless 219 (Sh219) has no dense 
shell . Only a diffuse neut ral layer surrounds the ionized 
gas l|Roger fc Leahvlll993|h though the sign of triggering 
is shown in t he adjacent molecular cloud (De harveng et alJ 
l2003all2006ll . We show that these observational properties 
are successfully explained by our model, where the HII 
region and PDR expand in a radially stratified molecular 
cloud. 

The structure of this paper is as follows: In § [2 we 
briefly probe the expansion of the H II region and PDR in 
the radial density gradient with an analytic treatment. We 
focus on the time evolution of the column densities of the 
shell and envelope, which are key physical quantities for 
the shielding of FUV radiation. In § [31 we investigate the 
dynamical evolution and trapping of the IF and DFs in 
the molecular cloud using detailed numerical calculations. 
We apply our numerical modeling to Sh219 in § 0J We 
assign §0and|n|to discussions and conclusions. 
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2. Analytic Treatment 

In this paper, we study the simultaneous expansion of the 
HII region and surrounding PDR in an ambient medium 
presenting a radial density gradient. We adopt a density 
distribution including a central core and outer envelope, 



and mean that the mass ratio, M sw (t)/Mi(t) is pro- 
portional to p(R)/pi(t) oc R(3~2w)/2 

Since the SF appears 

at Ri ~ R w , we get 



M sw {t) 
Mi(t) 



R(t) 



-i (3-2w)/2 



R-u 



(6) 



n(r) 



for r < R c 
for r > R c 



(1) 



Using equations 0, and 0, the HII mass density is 
written as, 



where R c is the core radius, n c is the core density, and 
w(> 0) is the power-law index of the density profile in 
the outer envelope (FTB90). First, we consider the expan- 
sion of the H II region in the analytic approach, following 
FTB90. The expansion is subject to the balance between 
the recombination rate in the HII region and supply rate 
of UV photons by the star. Unless this balance is achieved, 
the IF quickly propagates in the molecular cloud as an R- 
type front. It is not until the IF runs out of the cloud that 
the ionized gas expands owing to the pressure gradient 
(radial "champagne flow"). Once the UV photon supply is 
balanced by the recombination, on the other hand, the IF 
decelerates and is "trapped" in the cloud. The IF changes 
from R-type to D-type, and the SF emerges in front of the 
IF. FTB90 have shown that this can occur when w < 1.5. 
At the moment that the ionization balance is achieved, 
the radius of the HII region is, 



(3-2uj)/2 



(7) 



In this paper, we focus on the structure of the outer PDR 
as well as the HII region. The PDR structure depends 
on the efficiency of the shielding of FUV radiation from 
the central star. Since the FUV radiation is attenuated by 
the dust absorption and self-shielding effect, the penetra- 
tion of the PDR is sensitive to the column density of the 
shell and outer envelope. First, let us evaluate the time 
evolution of the column density of the shell. Calculating 
M svf (t) - Mi(t) with equations 0,0, and 0, the mass 
of the shell at the time t is, 



M sh (t) 



in 
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Therefore, the column density of the shell is, 
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for w < 1.5, and 



R1.5 = -R c exp 



#7 



(3) 



for w — 1.5, where R st is the Stromgren radius correspond- 
ing to the core density. Note that R w becomes equal to R st 
with w = 0. The subsequent expansion is driven by the 
high pressure of the H II region. A shell forms around the 
HII region. A part of the swept-up mass flows into the 
HII region, and the other part remains in the shell. The 
total swept-up mass at the time t is written as, 



Equation shows an interesting feature of the time evo- 
lution of the shell column density. As the H II region ex- 
pands, the shell column density increases with w < 1, 
but decreases with w > 1 . The critical power-law index is 
w = 1, with which the shell column density converges to 
the constant value of n c R c /2. If the density gradient is as 
steep as w ~ 1.5, the column density of the shell is signifi- 
cantly reduced. Supposing the finite size of the cloud, the 
inner edge of the envelope is the IF and the outer edge 
is the end of the cloud. The column density of this outer 
envelope is, 



M sw (i) 



4tt 



-p c Rl 



R(t) 



Airr dr 



N ev (t) = 



47T 

3 — w 



P (R)R(t) 3 
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where R(t) is the position of the SF, and p(R) is the initial 
mass density at r = R given by equation . The mass of 
the ionized hydrogen in the H II region is given by 
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w — 1 



n(r) dr 
Rc 



R(t) 
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Mi(t) = yPii?i(i) 3 



(5) 



where Ri(t) is the position of the IF, and pi(t) is the mass 
density in the HII region at the time, t. Since the total 
recombination rate within the H II region always balances 
the UV-photon number luminosity, the HII density de- 
creases as pi cc i?; -3 ^ 2 . Presuming R(t) ~ Ri(t), equations 



for w 7^ 1, where i? ut is the radius of the outer limit 
of the cloud. Equation i|10|) shows that the column den- 
sity of the envelope is mainly determined by i? ou t with 
w < 1, but R(t) with w > 1. Therefore, if the density 
gradient is as steep as w > 1 , iV cv quickly decreases as the 
HII region expands. Equations and (|10f) predict that 
the time evolution of the outer PDR qualitatively changes 
across w ~ 1. The quantitative changes will be signifi- 
cant, when the initial column density of the core is about 
7V H ~ 4 — 8 x 10 21 cm~ 2 . In a molecular cloud stratified 
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with w < 1, the column density of the shell increases as 
the HII region expands, and the dust absorption becomes 
significant in the shell. The molecules are shielded from 
the FUV radiation by the dust, and will accumulate in 
the shell. This is similar to the evolution in the homoge- 
neous ambient medium studied in Papers I and II. With 
w > 1, on the other hand, the shell column density finally 
decreases as the H II region swells up. The dust absorption 
does not work in the shell, and the PDR extends far be- 
yond the shell. Furthermore, equation IjlOjl shows that the 
column density of the outer envelope quickly decreases as 
the HII region expands. Therefore, it becomes easier and 
easier for the FUV photons to penetrate and photodissoci- 
ate the outer molecular envelope. In the following sections, 
we verify our predictions, performing detailed numerical 
calculations. 

3. Numerical Calculations 

In this section, we present our numerical results calculated 
with our code developed in Papers I and II (see Paper II 
for detailed description) , with which we can solve the time 
evolution of the thermal and chemical structure of the 
outer PDR as well as of the dynamical evolution of the HII 
region. The UV/FUV radiative transfer and some chemi- 
cal rate equations are consistently solved with the hydro- 
dynamics. The dust absorption is included only outside 
the HII region. We adopt the ambient density distribution 
given by equation with n c = 10 5 cm -3 , and the excit- 
ing star of 40.9 M©. The corresponding stellar UV- and 
FUV-photon number luminositi es are 5ttv = 6. Ox 10 48 s" 1 
and Sftjv = 5.8 x 10 48 s" 1 ijDiaz-Miller et al. Ill998h . 
The Stromgren radius with the adopted n c and Suv is 
Rst ~ 2.5 x 10~ 2 pc, and we set the core radius as 
R c = 0.7i? st ~ 1.75 x 10" 2 pc. The IF and the DFs are 
initially set at r = 0.2 R c and r = 0.4 R c respectively. 

Figures ^ an d 13 present the gas-dynamical evolution 
with w — 0.8 and 1.2 respectively. The dynamical fea- 
tures are similar in both models. When the IF reaches the 
radius, R w , the SF appears in front of the IF. FTB90 have 
provided the approximate equation for the time evolution 
of the H II radius in this phase, 



R(i) — R v 



7-2w 



12 



t 



9 - 4u> t d yn 



4/(7-2m>) 



(11) 



where tdyn is defined as R w /C\\, and Cjj is the sound 
speed in the HII region (T ~ 10 4 K). Equation still 
provides a good approximate evolution for our numerical 
calculations. In both models, the density of the shell is 
10-100 times as dense as that just ahead the SF. The ex- 
pansion is faster with the steeper density gradient. The 
size of the H II region with w — 1.2 is about twice as large 
as that with w — 0.8 at the same time, t. The outward 
velocity within the H II region is larger with the steeper 
density gradient. 

Despite their similar dynamical evolution, the thermal 
and chemical structures of these model's shells and outer 




Fig. 1. Gas-dynamical evolution with the density power- 
law index, w = 0.8. In each panel, we present five snap- 
shots at t = 0.033, 0.1, 0.167, 0.23, and 0.3 Myr. 



envelopes are very different. The third panels of Figures Q 
and [21 show the temperature profile in each model. With 
the density gradient of w — 0.8 (Fig^J, the temperature 
in the outer envelope is kept at T ~ 10 — 30 K, except for 
just outside the shell. With the steeper density gradient 
of w = 1.2 (FigEJ, however, the outer envelope is heated 
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Fig. 2. Same as Fig^but for the density profile with w 
1.2. 



up to T ~ 100 - 300 K in the initial 0.1 Myr. Figure El 
shows the chemical structure of the shell and outer enve- 
lope with w — 0.8 at each time step. The DFs of both 
H2 and CO molecules are trapped in the shell by the time 
of t ~ 0.1 Myr. The molecular gas gradually accumulates 
in the outer region of the shell. The outer envelope sur- 
rounding the H II region remains molecular. In the model 
with w — 1.2, on the other hand, wc can sec the different 




0.15 0.2 
Time (Myr) 

Fig. 3. Upper panel : Time evolution of some front po- 
sitions in the model with w = 0.8. The solid, broken, 
and dot-solid lines represent the positions of the IF, DF 
of the H2 and CO molecule at each time step. We do 
not present the position of the SF, which almost traces 
the position of H2 DF. Lower panel : The time evolution 
of the column density of each region in the same model. 
The upper (lower) solid line indicates the position of the 
SF (IF). Two broken (dot-solid) lines show the position 
where xh 2 (xco) = 0.1 and 0.9. The shaded region de- 
notes the expected gravitationally unstable region, where 

t > 1/y/Up. 

chemical evolution in the upper panel of Figure While 
the H 2 DF remains just in front of the SF, the CO DF 
leaves the SF and quickly travels over the whole cloud. 
Almost all CO molecules within 10 pc are photodissoci- 
ated by the time of t ~ 0.2 Myr. These clear differences 
are due to the different evolution of the column densities 
of the shell and outer envelope, which was predicted in 
§0 In the model with w = 0.8, the column density of the 
shell increases as the HII region expands, as the lower 
panel of Figure and Figure [S] show. The shell column 
density attains Ngt ~ 5 x 10 21 cm -2 , which corresponds 
to Ay ~ 2.5, at t ~ 0.1 Myr. The significant dust ab- 
sorption in the shell efficiently blocks FUV photons from 
the central sta r. The FUV radiation field in Habing units 
llHabindll96S|) is G FU v ~ 10 4 at the IF, but only - 50 at 
the SF at t ~ 0.1 Myr. Consequently, the gas temperature 
of the outer envelope is as low as T ~ 10 — 30 K, owing to 
the weak FUV radiation field. Both H2 and CO molecules 
are protected against FUV photons by the dust absorp- 
tion and self-shielding effect, and DFs are easily trapped 
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0.05 0.1 0.15 0.2 0.25 0.3 
Time (Myr) 



Fig. 4. Same as Fig|3]but for w = 1.2. We plot the posi- 
tion of the SF with a dotted line in the upper panel. Note 
the differences in scale of the vertical axis from FigEI The 
dotted line in the lower panel represents the column den- 
sity at r ~ 10 pc. Almost all CO molecules within 10 pc 
are photodissociated by the time of t ~ 0.2 Myr. 



in the shell. After the DF is engulfed in the shell, the SF 
sweeps up the molecular gas in the envelope. The molecu- 
lar gas accumulates in the shell; this is basically the same 
time evolution as in a homogeneous ambient medium (see 
Papers I and II). 

With the steeper density gradient with w = 1.2, our 
model shows the different time evolution of the column 
density in each region (the lower panel of Fig 0] and Fig|3|). 
As expected with equation ©, the column density of the 
shell decreases as the HII region expands. Since the shell 
column density never exceeds 1.0 x 10 21 cm -3 , the dust ab- 
sorption in the shell is not sufficient to protect molecules. 
The self-shielding effect barely enables the gradual accu- 
mulation of H2 molecules in the shell, but the FUV ra- 
diation easily photodissociates all CO molecules in the 
shell. Furthermore, the column density of the outer enve- 
lope also quickly decreases following equation (|10|) . The 
column density becomes too low to shield the FUV radi- 
ation from the central star by the dust absorption in the 
envelope. For example, the column density at r ~ 10 pc 
is initially higher than 2.0 x 10 22 cm~ 3 , but lower than 
4.0 x 10 21 cm" 3 at t ~ 3 Myr (FiggJ. The DF of the 
CO molecule leaves the SF, and travels over the whole 
molecular cloud. 
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Fig. 5. Time evolution of the column density of the shell 
with the different density distributions. The three solid 
lines represent the time evolution expected by equations 
© and (nj with w = 0.8, 1.0 and 1.2 from top to bot- 
tom. The open circles denotes the shell column density 
calculated every 0.05 Myr in our numerical models. The 
horizontal dashed line indicates the half column density of 
the initial homogeneous core, n c R c /2 ~ 2.7 x 10 21 cm~ 2 . 

In the cloud with a much steeper density gradient, it 
becomes much easier for the DFs to escape from the dense 
central region of the cloud. As referred to in § [3 the H II 
region can be trapped in the cloud with a density stratifi- 
cation of w < 1.5 (FTB90). Therefore, only the PDR can 
broadly extend around the trapped H II region with the 
density profile of 1 < w < 1.5. With such a steep density 
gradient, and if the DFs reach the diffuse outskirts of the 
cloud, a density-bounded PDR should be observed around 
the photon-bounded H II region. 

4. Application to the Galactic H II Region ; 
Sharpless 219 

4.1. A Model with w = 1.5 

In this section, we focus on the Galactic H II region, 
Sharpless 219 (Sh219). Sh219 shows a spherical morphol- 
ogy around a single B0V star; its radius is about 2.2 pc. 
As mentioned above, Sh219 shows some properties dif- 
ferent from those of the "collect and collapse" HII re- 
gions. For example, the expanding shell around the H II re- 
gion has not detected by the molecular emission. Only the 
half-ring of the PAH em i ssion surrounds the ionized gas 
llDeharveng et al.1 l2003at iDeharveng. Zavagno fe Caplanl 
l2005h . Instead of a dense molecular shell, a diffuse atomic 
layer is widely distributed around the H II region, and the 
diff use far-IR emission has been detected in this region 
l|Roger fc Leahvl [l993) . The radial width of the neutral 
layer is 2 — 3 pc. The averaged number density of the 
layer is only ~ 9 cm~ 3 , which is as low as the density 
typical of the cold neutral medium, rather than that of 
the molecular cloud. The number density in the HII re- 
gion has been estimated as 55 cm" " 3 l|Roger fc Leahvll993[) 
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Fig. 6. Snapshots of the chemical structure of our model 
for Sh219. The thin dot-solid and broken line represent 
the initial molecular ratios, Xqo = "■co/^H.nuc^C and 
Xh 2 = 2nH 2 /^H,nuc, where nn,nuc is the number density 
of the hydrogen nuclei. The thick solid and broken line 
represent the ratios, Ajjii = /«H,nuc and Ah 2 at t = 
0.01 Myr. 

and 170 cm IIDeharveng et al . 2000). Therefore, the ion- 
ized gas is much denser than the surrounding neutral 
layer. These features sharply contrast with those of a "col- 
lect an d collapse" H II region, fo r example, Sharpless 104 
(Shl04: iDeharveng et aLll2003rj|) . Shl04 is also a pc-scale 
spherical H II region, but surrounded by a dense molecular 
shell. The estimated average density and molecular mass 
of the shell are about 6000 cm -3 and 6000 Af Q respec- 
tively. The number density in the H II region is ~ 80 cm~ 3 , 
which is much lower than that of the shell. We have shown 
that these observational properties of Shl04 are well ex- 
plained by our model, where the H II region swells in the 
homogeneous ambient medium of ~ 10 3 cm~ 3 (Paper I). 
In this subsection, we propose an alternative model for 
Sh219, where the HII region is expanding in a radially 
stratified molecular cloud. 

We adopt the initial density distribution given by 
equation 0) with w = 1.5. We adopt a 19 Af© star as 
the central star, and the corresponding UV and FUV lu- 
minosities aje_j5j^_^_5 ; 6_xl0^_s _1 and Sfuv = 1-6 x 
10 48 s" 1 l|Diaz-Miller et al. I Il99^) . The time evolution 
quantitatively changes with the different parameters of the 
density and radius of the core. In this section, we present 
a model with n c = 10 5 cm~ 3 and R c = 0.7 R st ~ 0.07 pc. 
The observed neutral layer around Sh219 is so diffuse that 
we take into account the fact that the FUV background 
radiation (FBR) penetrates into the cloud. The initial 
thermal and chemical structure of the cloud are set at 
the equilibrium determined with the local number den- 
sity and FUV radiation field. We adopt Gfuv = 1-0 in 
Habing units as the fiducial FBR. Note that the strength 
of the FBR does not affect the chemical/thermal structure 
outside the H II region when the IF reaches the observed 
radius of Sh219 (see below). The third panel of Figure 



{7\ shows that the gas temperature in the outer envelope 
is initially as high as 100-200 K, owing to the photoelec- 
tric heating by the FBR. The photodissociation by the 
FBR is significant in the outer diffuse region. More than 
90% of H 2 (CO) molecules is initially photodissociated at 
r > 3 pc (r > 0.3 pc) (see FigEJ). Figure [7| shows the gas- 
dynamical evolution of our model for Sh219. The basic 
evolution is similar to the models presented in § The IF 
and preceding SF expand, and reach the radius of Sh219, 
r ~ 2.2 pc at t ~ 0.1 Myr. At that time, the number den- 
sity in the HII region is ~ 50 cm -3 , and that of the outer 
layer is ~ 10 cm~ 3 at r ~ 4 pc. Therefore, the observed 
density structure of the Sh219 is naturally explained by 
our model. 

The molecular gas in the central dense region is quickly 
destroyed by the FUV radiation from the central star. H 2 
molecules in the inner 1 pc are photodissociated over the 
initial 0.01 Myr (see FigHJ). The upper panel of Figure |S1 
presents the time evolution of the chemical structure of the 
envelope. As this panel shows, H 2 molecules in the enve- 
lope completely disappear by the estimated age of Sh219, 
t ~ 0.1 Myr. Almost all H 2 molecules are dissociated by 
~ 0.04 Myr, and CO molecules in the inner region are 
destroyed just after the start of the calculation. The H II 
region is surrounded with the diffuse atomic layer, where 
the gas temperature is about 300 — 500 K. Our model sug- 
gests that the observed wide neutral layer around Sh219 
should not be photon-bounded, but density-bounded. This 
does not depend on the strength of the FBR. In Figure 
|H1 we also plot the propagation of the H 2 DF when the 
FBR is not included. Even without the FBR, FUV pho- 
tons from the central star easily penetrate into the molec- 
ular envelope, and arc sufficient for the photodissociation 
of almost all H 2 molecules in the envelope. This is because 
the column densities of the shell and envelope are too low 
to block FUV photons with the steep density gradient of 
w = 1.5 (see §EJ. 

Figure [3] shows the time evolution of the mass of the 
HII region and outer neutral layer. Assuming that the 
neutral layer is density-bounded at r ~ 4 — 5 pc, which is 
the observed outer radius of the H I layer, the calculated 
masses at t ~ 0.1 Myr, when the IF reaches the observed 
radius of Sh219, show good agreement with the observed 
values. The mass of the H I layer, Mhi, initially increases, 
because the H 2 molecules in the envelope are photodisso- 
ciated by the FUV photons from the central star. After 
all H 2 molecules are destroyed (t ~ 0.02 Myr), Mm de- 
creases, because the H II region gradually erodes the H I 
layer. Note that the mass of the shell is much smaller 
than that of the H II region and neutral layer. Most of the 
swept-up gas does not remain in the shell, but flows into 
the HII region. 

4.2. Timescale of Induced Star Formation 

Although Sh219 does not have a clear molecular shell, 
some observations have suggested that star formation is 
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actually triggered i n a molecular cloud present at the 
periphery of Sh219 foeharveng. Zavagno fe Caplanl feoOS: 
iDeharveng et all l2006|) . A voung stellar cluster is emb ed- 
ded in the molecular cloud, and elongated along the IF. 
The cluster includes several massive stars exciting a ultra- 
compact H II region. One possible triggering scenario is 
that the IF and the preceding SF enter the pre-existing 
molecular cloud, and the dense compressed layer forms 
only in the cloud. The fragmentation of the compressed 
layer occurs and triggers star formation. Here, we exam- 
ine this scenario, evaluating the timescale of this triggering 
process. 

The average number density of the molecular cloud is 



estim ated as n(H 2 ) = 8.0 x 10 3 cm 3 1 Deharveng et alJ 



2006). In our model, the expanding velocity of the IF and 
SF is 2 chii ~ 20 km/s (chii is the sound speed in the 
HII region). After the IF and SF hit the molecular cloud, 
the SF can still propagate at ~ several - 10 km/s within 
the photo-evaporating cloud l|Bertoldil ll989). If the cloud 
size is much smaller than the H II radius, the SF quickly 
transverses and compresses the cloud. The Mach number 
of the SF will be of the order of 10, and the density jump 
at the isothermal SF is, 



n-y 



(12) 



where ni and n 2 are number densities ahead/behind the 
SF. With the average cloud number density, the number 
density just behind the SF can be n(H 2 ) ~ 10 6 cm~ 3 . The 
timescale of the fragmentation is, 



8.4xl0 4 y/ " (H2) 



V io 6 



-1/2 



(13) 



and tf rg with n(H2) ~ 10 6 cm -3 is comparable with the 
estimated age of Sh219, ~ 0.1 Myr. Therefore, triggering 
is possible. This is also comparable with the typical age of 
the u l tra-compact HII regions fe.g.. lWood fc ChurchweJ 
Il989t IChurchwel]ll2002h . Note that if Sh219 is modeled 
in a homogeneous ambient medium of 30 — 50 cm -3 , the 
estimated age is much younger, only a few times 10 4 yr 
I Roger fc Leahvl ll9931. In this case, some special condi- 
tions should be satisfied for the rapid triggering scenario: 
The SF must enter the region much denser than average 
density, and subsequent star (cluster) formation must ad- 
vance very rapidly. 



5. Discussions 

5.1. Implication for the Feedback Effect : Negative or 
Positive ? 

In this subsection, we discuss the role of UV/FUV 
radiation for star forma tion in the molecular cloud . 
iRoger fc Dewdnevl l|l992|) and biaz-Miller et al. I l|l998|) 
have calculated the expansion of the H II region and PDR 
in a homogeneous ambient medium, solving the UV/FUV 
radiative transfer around a massive star. They have shown 



that a significant amount of the mol ecular gas is pho- 
todiss ociated by the FUV radiation. iDiaz- Miller et al.l 
l|l998|) have suggested that the star formation efficiency 
of the molecular cloud is significantly suppressed by the 
photodissociation (negative feedback). However, our up- 
to-date calculations, which solve the hydrodynamics as 
well as the radiative transfer, have shown another aspect 
of the feedback effect (Papers I and II). We have shown 
that the dense shell forms, and that most of the swept- 
up gas remains in the shell in a homogeneous medium. 
While the FUV radiation destroys the molecular mate- 
rial, FUV photons are finally shielded by the high col- 
umn density of the shell. The dense molecular shell forms, 
where triggering of the subsequent star formation occurs 
(positive feedback). Our quantitative analysis has shown 
that the positive feedback effect can dominate the nega- 
tive effect in some cases with the homogeneous medium 
l|Hosokawa fc Inutsukall2006bj) . 

Our current work has shown that the negative feedback 
effect is promoted only in the inhomogeneous medium. 
In the cloud with a steep density gradient of w > 1, 
the shielding in the shell and envelope becomes ineffi- 
cient as the HII region expands. The PDR can widely 
extend around the ionized gas. Consequently, the star for- 
mation efficiency will be reduced, though triggering still 
can occur in the adjacent molecular clouds, as shown in 
Sh219. Ultimately, the feedback effect by the UV/FUV 
radiation should be examined in the clumpy, turbulent 
medium. Some recent studies have calculated the expan- 
sion of the HII region in this realistic situation (e.g., 
IMellema et al1l2006l:lMac Low et al.ll2006|) . and suggested 
that t he clumpy structu re may also diminish the positive 
effect l|Dale et al.ll2005T) . Although these efforts have not 
included the outer PDR, similar approaches solving the 
FUV radiative transfer will make it possible to clarify the 
role of the FUV radiation. For this purpose, however, the 
swept-up shell should be resolved with a sufficient number 
of grids. The column density of the shell often dominates 
the total column density of the system, which is important 
for the penetration of FUV photons. 



6. Conclusions 

In this paper, we have studied the time evolution of the 
HII region and surrounding PDR in the radially stratified 
molecular cloud. We have examined the efficiency of the 
trapping of FUV photons in clouds with different density 
profiles represented as n(r) cx r~ w . We have focused on 
the expansion with w < 1.5. In the molecular cloud with 
w > 1.5, neither the IF nor DFs are trapped, and quickly 
travel over the whole cloud. 

The key physical quantity for the trapping of the FUV 
radiation is the column density, because the dust absorp- 
tion is the primary shielding agent. First, we have analyt- 
ically shown that the time evolution of the column densi- 
ties of the shell and outer envelope qualitatively switches 
across w = 1. The column density of the shell increases 
as the HII region expands in the cloud with w < 1. With 
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w > 1, however, the shell column density finally decreases, 
and the column density of the outer envelope quickly de- 
creases as the HII region expands. The quantitative dif- 
ference across w = 1 is significant when the initial column 
density of the core is ~ several x 10 21 cm~ 2 . 

Next, we have verified the analytic consideration using 
the numerical calculations. The chemical/thermal struc- 
ture outside the HII region sharply changes with w. In 
the cloud with a steep density gradient of w > 1, the PDR 
can extend broadly around the HII region. The stellar 
FUV radiation heats up the envelope to several x 100 K 
via the photoelectric heating. The trapping of the DFs 
is sometimes selective. Only CO DF can travel over the 
whole cloud, while H 2 DF is trapped in the shell. In the 
steeper density gradient, the trapping of the DFs grows 
less efficient. If the density gradient is as steep asw~ 1.5, 
both H 2 and CO DFs quickly propagate in the cloud, and 
almost all molecules in the envelope are finally photodis- 
sociated. This is contrasted with the expansion in the ho- 
mogeneous medium and the gradual density gradient of 
w < 1. In these cases, DFs are finally engulfed by the 
shell, and the molecular gas gradually accumulates in the 
shell (Papers I and II). 

Finally, we have applied this case (w — 1.5) to 
the Galactic HII region, Sh219, whose observational 
properties are very different from the "collect and col- 
lapse" candidates. In Sh219, the ionized gas is not sur- 
rounded by a dense molecular shell, but by a diffuse neu- 
tral layer. The H II region is denser than the surrounding 
neutral layer. Our model naturally explains these charac- 
ters of Sh219. The calculated size, density and mass of 
the H II region and neutral layer are in good agreement 
with the observed values at t ~ 0.1 Myr. We suggest that 
a density-bounded PDR surrounds the photon-bounded 
HII region in Sh219. 
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Fig. 7. Gas-dynamical evolution of our model for Sh219. 
The adopted power-law index of the density profile is w — 
1.5. We present five snapshots at the time of t = 0.015, 
0.045, 0.075, 0.105, and 0.135 Myr in each panel. The ver- 
tical dot-solid line indicates the observed radius of Sh219, 
- 2.2 pc. 
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Fig. 8. Upper panel: Time evolution of the positions of 
some fronts in our model for Sh219. The solid and dot- 
ted line represent the positions of IF and SF. The broken 
contours denote the positions where the H2 molecular per- 
centages are 5%, 25%, 50%, and 75% from up to down. 
The thin broken line shows the position of the H2 DF with- 
out the FUV background radiation. The horizontal dashed 
line indicates the observed radius of Sh219, ~ 2.2 pc. The 
HII radius reaches the observed radius at t ~ 0.105 Myr, 
which is indicated with vertical dashed line. Lower panel: 
Time evolution of the mass in each region. The solid and 
dotted lines represent the mass of the ionized hydrogen 
in the H II region, and neutral hydrogen in the shell re- 
spectively. Two broken lines show the mass of the diffuse 
neutral layer in r < 4 and 5 pc. Asterisks indicate the ob- 
served mass of the HII r egion (34 Mr.-,) and sur rounding 
neutral layer (86 M Q ) bv lRoger fe Leahvl l)l993(l . 



